Energy dependence of elliptic flow from heavy-ion collision models 
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We have compared the experimental data on charged particle elliptic flow parameter (V2) in 
Au+Au collisions at midrapidity for -^/snn = 9.2, 19.6, 62.4 and 200 GeV with results from various 
models in heavy-ion collisions like UrQMD, AMPT, and HIJING. We observe that the average V2 
((1*2)) from the transport model UrQMD agrees well with the measurements at ^/snn = 9.2 GeV 
but increasingly falls short of the experimental (v?) values as the beam energy increases. The 
difference in (V2) being of the order of 60% at snn = 200 GeV. The (v?) results from HIJING 
is consistent with zero, while those from AMPT with default settings, a model based on HIJING 
with additional initial and final state rescattering effects included, gives a {V2) value of about 4% 
for all the beam energies studied. This is in contrast to increase in {1J2) with beam energy for the 
experimental data. A different version of the AMPT model, which includes partonic effects and 
quark coalescence as a mechanism of hadronization, gives higher values of (t^) among the models 
studied and is in agreement with the measured (V2) values at ^/snn = 200 GeV. These studies show 
that the experimental (V2) has substantial contribution from partonic interactions at ^/snn = 200 
GeV whose magnitude reduces with decrease in beam energy. We also compare the available data 
on the transverse momentum and pseudorapidity dependence of V2 to those from the above models. 
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Elliptic flow (1*2) measured in heavy-ion collisions 
are believed to arise because of the pressure gradient 
developed when two nuclei collides at non-zero im- 
pact parameters followed by subsequent interactions 
among the constituents {!]. Within a hydrodynami- 
cal framework, the vi has been shown to be sensitive 
to the equation of state of the system formed in the 
collisions @. Recent data on v 2 for baryons and 
mesons when measured as a function of tot — rrih, 
where tot (= \J P% + m \) is the transverse mass, 
Pt is the transverse momentum and to/j is the mass 
of the hadron, show a unique scaling at ^/snn = 
200 GeV. When v 2 and tot — win are scaled by the 
number of constituent quarks for a hadron, the v 2 
values follow an universal scaling for all the mea- 
sured hadrons This observation, refered to as 
the number of constituent quark scaling, has been 
interpreted as the collectivity being developed at 
the partonic stage of the evolution of the system 
in heavy-ion collisions [1, Further it has been 
shown that the pseudorapidity (77) dependence of v 2 
for charged hadrons shows a longitudinal scaling as 
observed for the multiplicity distributions in these 
collisions Is 01. 
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FIG. 1: (Color online) Average elliptic flow ((1*2)) as 
a function of beam energy. The results are shown for 
charged particles from RHIC experiments of STAR [l2T ]. 
PHENIX fp] and PHOBOS @|, SPS experiments of 
CERES [llfand NA49 (charged pions) [3], AGS experi- 
ments of E877 [H and E895 (proton) and FOPI 



The elliptic flow parameter is defined as the 2 nd 
Fourier coefficient v 2 of the particle distributions in 
emission azimuthal angle (<f>) with respect to the re- 
action plane angle (\f r ) [8|, and can be written as 



dN 



oc 1 + 2v 2 cos(2(0- *)). 



(1) 



For a given rapidity window the second coefficient is 

2 2 

« 2 = (cos(2(0-*))) = (%^%, (2) 

Px+Py 

where p x and p y are the x and y component of the 
particle momenta. 

The Relativistic Heavy Ion Collider (RHIC) has 
undertaken a beam energy scan program to look 
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for change in observations of various measurements 
as a function of beam energy to establish the par- 
tonic phase at higher energy collisions |9J. In such 
a program, the energy dependence of elliptic flow 
will play a crucial role. The compilation of available 
average charged particle (ua) as a function of beam 
energy [To] ] is shown in Fig. [TJ A non-monotonic 
dependence of (V2) vs. ^/snn is observed. 

At lower energies the negative (V2) is attributed 
to out-of-plane squeeze-out phenomena [13, U3 ■ I n 
this case the elliptical shape of the particle trans- 
verse momentum distribution at midrapidity is elon- 
gated in the direction perpendicular to the reaction 
plane and interpreted as due to shadowing by specta- 
tor nucleons. At high energies, the longitudinal size 
of the Lorentz contracted nuclei becomes negligible 
compared to its transverse size. This decreases the 
crossing time scales of the two nuclei. The shadow- 
ing effect goes away and elliptic flow fully develops 
in-plane, leading to a positive value of (1)2)- In this 
work, we concentrate on understanding this posi- 
tive (1*2) values at high energies measured at RHIC 
by comparing to available models of heavy-ion colli- 
sions. 

Various observables are compared to theoretical 
calculations to understand the physical mechanism 
behind the measurements. Some of the freque ntly 
used models in heavy- ion collisions are UrQMD [l9[ , 
AMPT [H and HIJING [Uj]. HIJING (or Heavy 
Ion Jet Interaction Generator) is an event gener- 
ator for heavy-ion collisions. It is a perturbative 
QCD inspired model which produces multiple mini- 
jet partons, these later get transformed into string 
configurations and then fragment to hadrons. The 
fragmentation is based on the Lund jet fragmenta- 
tion model [22j]. A parameterized parton distribu- 
tion function inside a nucleus is used to take into 
account parton shadowing. Such a model does not 
have the mechanism to generate elliptic flow, how- 
ever it would be interesting to know how much of 
the correlations among hadrons from the minijets 
contribute to the (v 2 )- 

The AMPT (A Multi Phase Transport model) 
uses the same initial conditions as in HIJING. How- 
ever the minijet partons are made to undergo scat- 
tering before they are allowed to fragment into 
hadrons. The string melting (SM) version of the 
AMPT model (labeled here as AMPT-SM) is based 
on the idea that for energy densities beyond a criti- 
cal value of ~ 1 GeV/fm 3 , it is difficult to visualize 
the coexistence of strings (or hadrons) and partons. 
Hence the need to melt the strings to partons. This 
is done by converting the mesons to a quark and 
anti-quark pair, baryons to three quarks etc. The 
scattering of the quarks are based on parton cas- 
cade [2(J. Once the interactions stop, the partons 
then hadronizes through the mechanism of parton 
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FIG. 2: (V2) for charged particles at midrapidity for 
minimum bias collisions at t/swn = 9-2, 19.6, 62.4 and 
200 GeV [3, E3] compared to corresponding results from 
AMPT and UrQMD model calculations. 



coalescence. The interactions between the minijet 
partons in AMPT model and those between partons 
in the AMPT-SM model could give rise to substan- 
tial (i>2)- Agreement between the data and the re- 
sults from AMPT-SM would indicate the contribu- 
tion of partonic interactions to the observed (v 2 ). 

The UrQMD (Ultra relativistic Quantum Molecu- 
lar Dynamics) model is based on a microscopic trans- 
port theory where the phase space description of the 
reactions are important. It allows for the propaga- 
tion of all hadrons on classical trajectories in combi- 
nation with stochastic binary scattering, color string 
formation and resonance decay. It incorporates 
baryon-baryon, meson-baryon and meson-meson in- 
teractions, the collisional term includes more than 
50 baryon species and 45 meson species. The com- 
parison of the data on (z^) to those obtained from 
UrQMD model will tell about the contribution from 
the hadronic phase. 

In this paper we compare the (112) for charged par- 
ticles at midrapidity for ^/s NN = 9.2, 19.6, 62.4 and 
200 GeV with those from UrQMD (ver.2.3 ), AMPT 
(ver. 1.11) and HIJING (ver. 1.35) models. The 
parton-parton interaction cross section in the string 
melting version of the AMPT is taken to be 10 mb. 
We also discuss the pr and 77 dependences of (1)2). 

Figure [5] shows the (1*2) for charged particles at 
midrapidity for various ^/snn for minimum bias (0- 
80%) collisions ij. The results for ^/s™ = 9.2 GeV 
are for minimum bias 0-60% collisions j9|. The (V2) 
value increases linearly from about 3% at 9.2 GeV 
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to about 5% at 200 GeV. The experimental data are 
compared to (^2) calculated from UrQMD, AMPT 
and AMPT-SM with default settings. The centrality 
selection is same for data and the models. In con- 
trast to observations from the data, the model (v 2 ) 
values either remain constant or decreases slightly 
with increasing ^/snn- The (v 2 ) value from UrQMD 
at 9.2 GeV and those from AMPT-SM at 200 GeV 
are in good agreement with the data. The (^2) values 
from AMPT lie intermediate to those from UrQMD 
and AMPT-SM. If we assume the (v 2 ) values from 
UrQMD to be the contribution from hadronic phase, 
then this contribution (v^ r ® MD /v 2 ata ) varies from 
100% to about 40% of the measured (v 2 ) as beam 
energy increases from 9.2 GeV to 200 GeV. The 
higher values of (v 2 ) in data indicate the possible 
contribution that can come in such transport mod- 
els due to inclusion of initial/final state scattering 
effects and/or due to partonic interactions. Com- 
parison with AMPT-SM reflects that at 62.4 and 
200 GeV, the (v 2 ) has contributions from partonic 
interactions. We have estimated the (1*2) values from 
HIJING and found them to be consistent with zero. 
The (v 2 ) values obtained from HIJING are -0.00006 
± 0.0004, 0.0001 ± 0.0005 and -0.002 ± 0.0005 at 
200, 62.4 and 19.6 GeV respectively. 

We now present the differential v 2 by comparing 
the px and r\ dependence from data with those from 
the models. Figure |3] shows the comparison between 
experimental data on vi vs. px at midrapidity for 
Au+Au collisions and models at ^/snn = 200 and 
9.2 GeV. The data for 200 GeV are for minimum 
bias 0-80% collisions while those for 9.2 GeV are for 
minimum bias 0-60% collisions. Reasonably good 
agreement between data and AMPT-SM is seen at 
200 GeV, whereas the AMPT-SM model results are 
higher at 9.2 GeV. This is consistent with the re- 
sults of (v 2 ) vs. y / SNN as shown in Fig. [2j Mea- 
surements exist for (i^) vs. px at midrapidity for 
Au+Au collisions at Jsnn = 62.4 GeV, the com- 
parison with models is consistent with the trends as 
revealed by the energy dependence shown in Fig. [51 
hence not presented here. The results from UrQMD 
and AMPT show similar trend of v 2 increasing with 
Pt but the values are lower compared to data at 200 
GeV. These models however are in a reasonably good 
agreement with the data at 9.2 GeV. For UrQMD, 
this is consistent with the results shown in Fig. O 
The AMPT results on (v 2 ) at 9.2 GeV was found 
to be higher than the data (Fig. [5]). The (v 2 ) folds 
the measured v 2 vs. px with the px distribution of 
charged particles. We found the (px) of charged par- 
ticles from AMPT is larger than those from UrQMD, 
resulting in the (v 2 ) being larger (as seen in Fig. [5]). 

FigureHJshows the comparison of the (^2) vs. 77 for 
charged particles (integrated over px) from Au+Au 
0-40% collisions at = 19.6 and 200 GeV [1. 
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FIG. 3: (Color online) v% as a function of pt at midra- 
pidity for Au+Au collisions at y/srra = 200 12] and 9.2 
GeV [IJ. The experimental results are compared to cor- 
responding vi values from UrQMD and AMPT models. 
The errors on the data points are statistical. 



The measurements show a decrease in (v 2 ) as the 
I 77 I increases. The general trend of (v 2 ) vs. 77 at 
62.4 GeV is also similar although the values lie inter- 
mediate to those from 19.6 and 200 GeV collisions. 
All the models also qualitatively shows similar de- 
creasing trend of (v 2 ) with increase in | 77 |, however 
they differ from the data in the widths of the (v 2 ) 
distribution in 77 and the values of (v 2 ). AMPT-SM 
model agrees with the data at midrapidity for y / snn 
= 200 GeV and is higher than the data at midra- 
pidity for ^/s NN = 19.6 GeV (consistent with results 
in Fig. However, for the higher rapidities the 
model fails quantitatively and has a larger width of 
the (^2) vs. 77 distribution compared to data. The 
AMPT and UrQMD models fail to explain quanti- 
tatively the dependence of (v 2 ) on 77 for ^/snn = 
200 GeV. At Js^ = 19.6 GeV the AMPT model 
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FIG. 4: (Color online) {V2) of charged particles as a func- 
tion of T) for Au+Au 0-40% collisions at ^/sNN = 19.6 and 
200 GeV [(J. The experimental results are compared to 
corresponding V2 values from UrQMD and AMPT mod- 
els. The errors on the data points are statistical. 



is closer to the data values for all the 77 range. The 
statistical errors on the data points are quite large. 

In summary, we have presented a compilation of 
the available data for elliptic flow parameter, V2, of 
charged particles at RHIC as a function ^snn, Pt 
and 77. These results are compared to correspond- 
ing model calculations from AMPT (default and 
string melting versions), UrQMD and HIJING. The 
(1)2) values at midrapidity increases with increase 
in »/snn- The (1*2) values from HIJING are con- 
sistent with zero at midrapidity for ^/snn = 19.6, 
62.4, and 200 GeV. This suggests that correlations 



among hadrons from minijets contribute negligibly 
to the actual (i^). The (1*2) results for the other 
models studied do not show a strong increase with 
V^nn- The AMPT-SM model agrees quite well with 
the (1*2) data at 200 GeV and gives higher values of 
(1)2) compared to data for the other beam energies 
studied. The results from UrQMD match with the 
(V2) data at 9.2 GeV and gives lower values of (^2) 
compared to data for the other beam energies stud- 
ied. The calculated (^2) values from AMPT model 
lie in between. Comparison with AMPT shows the 
growing importance of initial and final state scat- 
tering to (V2) with increase in ^/snn- Considering 
the (7J2) values from UrQMD to be the contribu- 
tion from hadronic matter, the higher values from 
the AMPT-SM model reflects additional contribu- 
tion due to partonic inter actions. The difference be- 
tween UrQMD and data (^2) values at midrapid- 
ity is about 60% at 200 GeV. This together with 
agreement of the experimental measurements with 
AMPT-SM results at 200 GeV clearly indicates sub- 
stantial contribution to the (1)2) measured at RHIC 
from partonic contributions. The agreement of mea- 
sured (^2) values with those from UrQMD at midra- 
pidity for 9.2 GeV suggests that the experimental 
results can be understood within the frame work of 
a hadronic model. 

We also presented a comparison of the measured 
i>2 vs. pt and rj with the various models. The mod- 
els qualitatively reproduce the trends of V2 values 
increasing with increase in pt and decreasing with 
increase in | r\ |. However, they lack in quantitative 
agreement at 200 GeV for the full rj range. Within 
the large statistical errors of the measured i>2 vs. r\ 
at 19.6 GeV, we observed the AMPT model gives 
the best description of the data. The RHIC beam 
energy scan program will make available the i>2 mea- 
surements at energies of 7.7, 11.5, 18, 27 and 39 GeV. 
One can then further systematically study the con- 
tribution of hadronic matter, initial state scattering 
and partonic matter to the V2 measurements through 
comparisons with the above models discussed. 
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